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Abstract 

Relativistic complex Burgers-Schrodinger and Nonliner Schrodinger 
equations are constructed. In the non-relativistic limit they reduce to 
the standard Burgers and NLS equations respectively and are integrable 
at any order of relativistic corrections. 

1 General Burgers-Schrodinger Hierarchy 

The relativistic linear Schrodinger equation has been discussed at the early years 
of quantum mechanics but was dismissed promptly by the Klein-Gordon and 
the Dirac equations. Recently, relativistic versions of the Schrodinger equation 
have been considered in the study of relativistic quark-antiquark bound states 
[1], and gravitational collapse of a boson star [2]. A nonlinear version of the 
model has appeared in the form of semi- relativistic Hartree-Fock equation [3]. 
But none of those models is known to be integrable. In the present paper we 
construct an integrable relativistic nonlinear Schrodinger equation, preserving 
integrability at any order of 1/c approximation. 

We start from the Schrodinger equation in 1 -I- 1 dimensions 

in^ = nPi)-^ (1) 

for a free particle with classical dispersion of the general analytic form E = 
E{p). Here Pq = ^nd Pi — —iti^ are operators of the time and space 
translations respectively, commuting with the Schrodinger operator S = — 
TC{Pi): [Pf^,S] — 0, ^ — 0,1. The general boost operator, defined as K = 
X — tTL^Pi), is also commuting with 5, [i^, S'] = 0. Commuting it with space 
and time translations we have the algebra of symmetry operators 

[Po,A] = o, [P(,,K]^~hn'{Pi), [Pi,K] = -ih. (2) 
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Then, if is a solution of ([T]) and W is an operator from this algebra, so that 
[W, S*] = 0, then S"^ is also solution of (P). 

For a given classical dispersion E = E{p), Eq = E{0) we define S-polynomials 
Hn {x,t) by the generating function 



^j-ipx-{E{p)-Eo)t) ^ V f - 1 ?- 



^h;-)(,„. (3) 

71=0 ^ ' 



It is equivalent to 

H\^\x,t) = ^-\{n{^-^dld.)-E,)t^n 
SO that Hn is a solution of 

th-^^Hi^\x,t) = {n- E,)Hi^\x,t) (5) 

( E) 

with the initial value Hn (x, 0) = a:". From commutativity [-S, X] — 0, time 
evolution of the operator K satisfies 

dK 

ih^ = [n,K] (6) 

and has the form K{t) = e"-E'"*ir(0) e^^* = e"'^'"*^ e*"^*. Then as follows, 
operator K generates the infinite hierarchy of polynomials according to 

KHnix,t) = ife-^("-^°)*a;" = e-i(^--f^«)*a;"+i = Hn+i{x,t) (7) 
1.1 Non-relativistic Schrodinger equation 

The non-relativistic dispersion E(p) = /2m implies the Hamiltonian operator 
and the Gallilean boost operator 



^ = - — 7^2 (8) 



K = x + it— — . (9) 

m ox 



From the generating function 

ib^-i4*) = ^(^) P-Hi^\x,t) (10) 



we have the Schrodinger polynomials 

i7^f)(a;,t) =e^**^^x". (11) 

If hIi^^^ {x, t) = exp{t-^)x'^ is the Kampe de Feriet polynomial, then iJn'^'' (x, t) - 
H^^{x, -^t) or in terms of the Hermit polynomial 



H^''Hx,t)= [-i^t] g„ I = ■ (12) 



2 



1.2 Semi-relativistic Schrodinger equation 



The relativistic dispersion E{p) — ^ vn?-d^ + (p-'p^ implies the Hamihonian 



and the semi-relativistic boost operator 



^ = 3337^2 (13) 



ih 



dx 



K^x + —t , (14) 

TO /-I 



In the non-relativistic limit, when c — > 00, it reduces to the Galilean boost ([9|). 
The generating function in the form of relativistic plane wave 



then gives semi-relativistic polynomials 

i) = . (16) 

In the non-relativistic limit H^^^ — > H^. The first three polynomials coincide 
exactly with the Schrodinger polynomials Hi^^ {x,t) = x, H2^^^\x,t) = x^ + 
i—t, iJo"^^"^' — x^ + i—ixt, while starting from the fourth one 

ilf (x, i) - + i-&xH -^3t^ + i^3t (17) 

TO TO^ TO-'C'^ 

we have relativistic corrections of order For complex valued space coor- 

dinate X, as it appears in 2+1 dimensional Chern-Simons theory [4], zeros of 
these polynomials describe a motion of point vortices in the plane. Equations 
of motion for iV vortices are (k = 1,...,N): 

1.3 Relativistic Burgers-Schrodinger Equation 

Using Schrodinger's log ^ transform 5 , '5 = e'"*, and identity 

the Schrodinger equation ([T]) can be rewritten in the form 
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For complex function = e's^^ = e''^"'"'^'^ we introduce a new complex function 
V = -i^-^ln^ = ^-^Py "^ith dimension of velocity. Then F = S - ihR, is 
the complex potential, real and imaginary part of which are the classical and 



quantum velocities V = Vc + iV„ = —S^ 



-R^. Hence (|20l) becomes of the 



complex Hamilton- Jacobi form (quantum Hamilton-Jacobi equation) 



F^]-1 = 



(21) 



In the classical (dispersionless) limit h ^ 0, the quantum velocity Vq vanishes 
and complex potential F reduces to the real velocity potential S, playing the 
role of Hamilton's principal function. In this case pT|) becomes the classical 
Hamilton-Jacobi equation ^ + 'H (ff ) — 0. By differentiation of ([20]) we have 
equation for the complex velocity 



ih 



dV_ 

dt 



h d 
m dx 



n -ih 



d_ 

dx 



mV\-l 



(22) 



This equation is the Madelung fluid type representation of the Schrodinger 
equation JT]). In the classical limit it gives the Newton equation rn^^ = 
~'§x [^(™^c)] or in the hydrodynamic type form 







(23) 



which is just differentiation of the classical Hamilton-Jacobi equation. Equation 
has implicit general solution as Vc(x,t) = f{x — T-C'{mVc)t), where / is an 
arbitrary function, and it develops shock at a critical time when derivative (Vc)x 
is blowing up. 



1.3.1 Non-relativistic Burgers-Schrodinger equation 

In this case the Schrodinger equation ([1]) with non-relativistic Hamiltonian 
is equivalent to the nonlinear equation for complex velocity V 



h^ d^v 

2m dx^ 



OX 



(24) 



which we call the Burgers-Schrodinger equation. In terms of the real and imag- 



inary parts it gives the Madelung fluid with density p 



and velocity Vc- In 



the classical limit it reduces to the one real equation for classical velocity Vc, 
namely the ordinary dispersionless Burgers equation. 



1.3.2 Semi-relativistic Burgers-Schrodinger equation 

For Hamiltonian ([13]), the " Burgerization" procedure described above gives the 
semi-relativistic Burgers-Schrodinger equation 



ldV_ d_ 
c dt dx 



1 



'ih— 
ox 



+ mV] • 1 



= 



(25) 
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In the non-relativistic limit it reduces to (23]) with rclativistic corrections in the 
lowest order as 



ih 



dV 
'dt 



ihv 



1 



2m dx^ 



dx 8m?c} 



[-h^V^xxx] (26) 



1 



+ — ^ [-imti^{lOV^V^^ + 41/T4,.) + rr?h'{l2VV^ + 6V^V,,) + Aim'^hV^V^] (27) 

In the classical (dispersionless) limit (j25p becomes equation of the hydrodynamic 
type 



iVc)t 



14 



(14). = 



(28) 



In the non-relativistic limit it reduces to the dispersionless Burgers equation 
with the lowest relativistic correction 



iVa)t + 14(K)x - ^14'(14):. = 



The general implicit solution of ([^5]) is 
Vcix,t) ^ f (^x ~ 
and it develops shock at a finite time. 



Vet 



(29) 



(30) 



1.4 Backlund Transformation for Burgers-Schrodinger equa- 
tion 

By the boost transformation of Section 1, from a given solution "ifi of the 
Schrodinger equation ([T|) we can generate another solution as 



Using identity 



*2 = K'i'i 



d 



x-tn' i -ih 



d_ 

dx 



dx 



^-^G -ih— W = G -ih— + mV\ - l 



d 



dx 



(31) 



(32) 



for complex velocities 14 = ^*:^ln.^a, (a = 1,2), we obtain the Backlund 
transformation 



h d 
m ox 



x-tH' { -ih-^ +mV] - l 
dx 



(33) 



For the non-relativistic quantum mechanics ([H]) it gives complex Backlund 
transformation 

V, = V,..1LI^ (34) 

m X — V^t 
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for the Burgers-Schrodinger equation (|24p . 

For the semi-relativistic quantum mechanics (|13[) we have the Backlund 
transformation of the form 

y2 = T/i-iAj-in|a;_ ^ v^t] (35) 

It is worth to note that in the classical limit ft 0, Vc the above Backlund 
transformations reduce to the trivial identity Vci — Vc2- 



2 Integrable General NLS Hierarchy 

In previous sections we studied the so called C-integrable relativistic Burgers- 
Schrodinger equation. Now using the AKNS hierarchy for the NLS equation we 
are going to construct relativistic NLS. 



2.1 NLS hierarchy 

We consider the Zakharov-Shabat linear problem 

d f vi \ _ f —^p —K^-il) \ f Vl \ _ ( Vl 



dx \ V2 J V V iP J \ '"'^ J V ^2 y 

for the space evolution, and the generalized AKNS problem [6] 

for the time evolution, where for the real A{x^t,p) and complex C{x,t,p) 
functions, determined by the zero-curvature condition, we substitute An = 
En=o^^"-' (-f)"' C!n = E,T=oC^"^ (-!)"• It gives the evolution equation 
dt^iP = d^C^°^ + 2iA^"^ and C^^) = 0, A^'^l = gn = const.. We fix 
this constant so that qn = (—2)^^^. Then we have the recurrence rela- 
tions C(") = ^9^C("+i) + ^("+1)7/;, (9^A(") = iK^j-^-Cn)^ _ c-l")^)^ where 
n = 0, 1, 2, N — 1. Integrating the last equation one has 

= -iK^ f (V^C(") - VC^")) (38) 



Substituting psp into recursion formula wc find 



where TZ is the matrix integro-differential operator - the recursion operator of 
the NLS hierarchy ^6^ - 

+ 2kV r V' -2AtV ^ 
7^ = ifT3 I I (40) 
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and (T3 - the Pauli matrix. Then we get 



-3( t) =7^-^ i ] (41) 



where In, N — 1,2,3,... is an infinite time hierarchy. In the hnear approx- 
imation, when K = 0, the recursion operator is just the momentum operator 
T^o = i'^3'§^ ^■iid the NLS hierarchy (|¥T|) becomes the hnear Schrodinger hierar- 
chy 

iil,^^ = (42) 

from Section 1. The Madelung representation for this hierarchy, produced by the 
complex Cole-Hopf transformation, is given by the complex Burgers hierarchy 

HI- 

Every equation of hierarchy (I4ip is integrable. The hnear problem for the 
A''-th equation is given by the Zakharov-Shabat problem (l36|) for the space part 
and 



dtM \ '"■2 J \ Cn -iAn J \ V2 J " \ V2 

for the time part. Coefficient functions Cn can be found conveniently as 

g ) =Ep-^-7^-^ ( ^ ) ^(p--+p--7^ + ... + 7^--) ( t 

k—1 

(44) 

To rewrite this expression in a compact form we introduce notation of the q- 
number operator 

l + q + q' + ... + q^-^ = [N]g (45) 

where g is a linear operator. Hence, with operator q = TZ/p we have following 
finite Laurent form in the spectral parameter p 



+ +■••+(-) ^[^Wp (46) 



Then we have shortly 



In a similar way 



g =."-MiVl«/,n) (47) 



•4» = - (/" *. - /' *) \N\ni, ( I ) (48) 

Equations , ([TT]) and give the time part of the linear problem (the Lax 
representation) for the N-th flow of NLS hierarchy (|4T|) in the q-calculus form. 
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2.2 General NLS hierarchy equation 

For time t determined by the formal series dt = '^'^^QE^dt^ where Ej^ are 
arbitrary constants, the general NLS hierarchy equation is 

za,^j^^^{Eo + E,n + ... + Er,n^ + ...) t ^ (49) 

Integrability of this equation is associated with the Zakharov-Shabat prob- 
lem ((551 and the time evolution 

Jo = £ E^Jo, = ( ) (50) 

where 



C 
C 



i:^»(g) = f:'^«''-iivw,(|) 

AT—r\ \ / AT—^ \ ^ / 



(51) 



N=0 ^ ' N=\ 

In the last equation we have used that for = 0, Cq = 0. Then we have 



N=0 N=0 ^ J J \ J 

The above equation (|49[) gives integrable nonlinear extension of linear Schrdinger 
equation with general analytic dispersion considered in Section 1. Let one con- 
siders the classical particle system with the energy-momentum relation E{jp) = 
Eq + Eip -\- E2P^ + .... Then the corresponding time-dependent Schrodinger 
wave equation is ([1]) where the Hamiltonian operator results from the standard 
substitution for momentum p —>■ —ifi-^ in the dispersion. Equation ([T]) together 
with its complex conjugate can be rewritten as 

-4(I)^-Hs)(l) 

The momentum operator here is just the recursion operator in the linear ap- 
proximation TZq = jfa^- Hence, (1531) is the linear Schrodinger equation with 
arbitrary analytic dispersion. The nonlinear integrable extension of this equa- 
tion appears as (HH), which corresponds to the replacement TZq ^ TZ, {h = 1), 
so that 



From this point of view, the standard substitution for classical momentum p —^ 
— ift^ or equivalcntly p — > —ifio's-^ = TZq for the equation in spinor form, gives 
quantization in the form of the linear Schrodinger equation. While substitution 
p ^ Ti- gives " nonlinear quantization" and the nonlinear Schrodinger hierarchy 
equation. 
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The related Lax representation for equation is given by ([?T|) . ([5^ . By 



definition of q-derivative Dq''^ f{() = '^^(gli for operator (7 = TZ/p, we have 
relation i*^/ C''^ = [^]tz/pP'^~^- Then equation ((5T|) can be rewritten as 



or using linearity of q-derivative and analytic dispersion form 

Due to above definition it gives simple formula 

C \ E{n) ~ E{p) ( V 



C J TZ-p V ^ 



Then for A we obtain 



lz..„x ..2^r.7 f\^\ E{n) - E{p) f 



(57) 



Equations ([57)1 . ([55]) give the Lax representation of the general integrable NLS 
hierarchy model (|54p in a simple and compact form. It is worth to note that 
special form of the dispersion E — E{p) is fixed by physical problem. In the next 
Section we will discuss the relativistic form of this dispersion and corresponding 
semi-relativistic NLS equation. 



3 Semi-relativistic NLS 

In Section 1.2 we have considered the relativistic dispersion relation E{p) = 
\J rn?c^ + p^c^. This may be used to construct a "semi-relativistic" Schrodinger 
equation with Hamiltonian (|13p . Then, combining two complex conjugate equa- 
tions together, we have 

We like to emphasize that if ■0 describes the relativistic particle forward in time 
and with positive energy, -0 corresponds to the backward time or to the negative 
energy. From this point of view equation (j59p is complete since includes both 
states. 

Following the general procedure described in previous Section one may pro- 
ceed further: by replacing the derivative operator TZq = *^3^ corresponding 
to linear momenta p with the full recursion operator TZ (|40p . one obtains an 
integrable relativistic nonlinear Schrodinger equation 
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where the square root operator has meaning of the formal power series so that 

V A V 2to2c2 8rn4c4 IGm^c^ J \ ^ 

(61) 

For the above relativistic dispersion and equation (|60p , we have the next hnear 
problem 



d 



Vl 



dx V ^^2 



d 



where 



-IP 








IP ) 


V V2 


-iA 




f"' 


C 


-iA ) 




TVc? - 







dt \ V2 



C J TZ-p V ^ 



(62) 
(63) 

(64) 



A = --v/m2c4+p2c2_iK2 / y V,-y V') (J 

(65) 

and the spectral parameter p has meaning of the classical momentum. The 
model (j60p . is an integrable nonlinear Schrodinger equation with relativistic 
dispersion: 

I 1 ^2~ 

^^t=mc\l--^—i: + F{i>) (66) 



, to2c2 (5a;2 
where the nonlinearity expanded in l/c2 is the infinite sum 



8m^c2 



It is interesting to note that if we expand also the dispersion part in l/c2, 
then at every order of l/c^ we get an integrable system. It means that we obtain 
integrable relativistic corrections to the NLS equation at any order. From the 
known relativistic integrable models like the Sine-Gordon or the Liouville equa- 
tions, neither one has this property. Finally we note that nonlinear relativistic 
equations considered in this paper are distinct from those obtained in [I]- [3] and 
references therein. They might be useful in analyzing relativistic corrections to 
solitons, Bose-Einstein condensates or other condensed matter systems with an 
efltective equation of relativistic form. 
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